Emergence of zoonotic sporotrichosis in Brazil: a genomic epidemiology study.

BACKGROUND: Zoonotic sporotrichosis is a neglected fungal disease, whereby outbreaks are primarily driven by Sporothrix brasiliensis and linked to cat-to-human transmission. To understand the emergence and spread of sporotrichosis in Brazil, the epicentre of the current epidemic in South America, we aimed to conduct whole-genome sequencing (WGS) to describe the genomic epidemiology. METHODS: In this genomic epidemiology study, we included Sporothrix spp isolates from sporotrichosis cases from Brazil, Colombia, and the USA. We conducted WGS using Illumina NovaSeq on isolates collected by three laboratories in Brazil from humans and cats with sporotrichosis between 2013 and 2022. All isolates that were confirmed to be Sporothrix genus by internal transcribed spacer or beta-tubulin PCR sequencing were included in this study. We downloaded eight Sporothrix genome sequences from the National Center for Biotechnology Information (six from Brazil, two from Colombia). Three Sporothrix spp genome sequences from the USA were generated by the US Centers for Disease Control and Prevention as part of this study. We did phylogenetic analyses and correlated geographical and temporal case distribution with genotypic features of Sporothrix spp isolates. FINDINGS: 72 Sporothrix spp isolates from 55 human and 17 animal sporotrichosis cases were included: 67 (93%) were from Brazil, two (3%) from Colombia, and three (4%) from the USA. Cases spanned from 1999 to 2022. Most (61 [85%]) isolates were S brasiliensis, and all were reported from Brazil. Ten (14%) were Sporothrix schenckii and were reported from Brazil, USA, and Colombia. For S schenckii isolates, two distinct clades were observed wherein isolates clustered by geography. For S brasiliensis isolates, five clades separated by more than 100 000 single-nucleotide polymorphisms were observed. Among the five S brasiliensis clades, clades A and C contained isolates from both human and cat cases, and clade A contained isolates from six different states in Brazil. Compared with S brasiliensis isolates, larger genetic diversity was observed among S schenckii isolates from animal and human cases within a clade. INTERPRETATION: Our results suggest that the ongoing epidemic driven by S brasiliensis in Brazil represents several, independent emergence events followed by animal-to-animal and animal-to human transmission within and between Brazilian states. These results describe how S brasiliensis can emerge and spread within a country. FUNDING: Coordenação de Aperfeiçoamento de Pessoal de Nível Superior, Brazil; the São Paulo Research Foundation; Productivity in Research fellowships by the National Council for Scientific and Technological Development, and Ministry of Science and Technology of Brazil.

Tracing histoplasmosis genomic epidemiology and species occurrence across the USA.

ABSTRACTHistoplasmosis is an endemic mycosis in North America frequently reported along the Ohio and Mississippi River Valleys, although autochthonous cases occur in non-endemic areas. In the United States, the disease is provoked by two genetically distinct clades of Histoplasma capsulatum sensu lato, Histoplasma mississippiense (Nam1) and H. ohiense (Nam2). To bridge the molecular epidemiological gap, we genotyped 93 Histoplasma isolates (62 novel genomes) including clinical, environmental, and veterinarian samples from a broader geographical range by whole-genome sequencing, followed by evolutionary and species niche modelling analyses. We show that histoplasmosis is caused by two major lineages, H. ohiense and H. mississippiense; with sporadic cases caused by H. suramericanum in California and Texas. While H. ohiense is prevalent in eastern states, H. mississipiense was found to be prevalent in the central and western portions of the United States, but also geographically overlapping in some areas suggesting that these species might co-occur. Species Niche Modelling revealed that H. ohiense thrives in places with warmer and drier conditions, while H. mississippiense is endemic to areas with cooler temperatures and more precipitation. In addition, we predicted multiple areas of secondary contact zones where the two species co-occur, potentially facilitating gene exchange and hybridization. This study provides the most comprehensive understanding of the genomic epidemiology of histoplasmosis in the USA and lays a blueprint for the study of invasive fungal diseases.

Bimodal distribution of azole susceptibility in Sporothrix brasiliensis isolates in Brazil.

Sporothrix brasiliensis is an emerging zoonotic fungal pathogen that can be difficult to treat. Antifungal susceptibility testing was performed on the mold phase of a convenience sample of 61 Sporothrix spp. isolates from human and cat sporotrichosis cases in Brazil using the Clinical and Laboratory Standards Institute standard M38. A bimodal distribution of azole susceptibility was observed with 50% (28/56) of S. brasiliensis isolates showing elevated itraconazole minimum inhibitory concentrations ≥16 µg/mL. Phylogenetic analysis found the in vitro resistant isolates were not clonal and were distributed across three different S. brasiliensis clades. Single nucleotide polymorphism (SNP) analysis was performed to identify potential mechanisms of in vitro resistance. Two of the 28 resistant isolates (MIC ≥16 mg/L) had a polymorphism in the cytochrome P450 gene, cyp51, corresponding to the well-known G448S substitution inducing azole resistance in Aspergillus fumigatus. SNPs corresponding to other known mechanisms of azole resistance were not identified in the remaining 26 in vitro resistant isolates.

Finding Candida auris in public metagenomic repositories.

Candida auris is a newly emerged multidrug-resistant fungus capable of causing invasive infections with high mortality. Despite intense efforts to understand how this pathogen rapidly emerged and spread worldwide, its environmental reservoirs are poorly understood. Here, we present a collaborative effort between the U.S. Centers for Disease Control and Prevention, the National Center for Biotechnology Information, and GridRepublic (a volunteer computing platform) to identify C. auris sequences in publicly available metagenomic datasets. We developed the MetaNISH pipeline that uses SRPRISM to align sequences to a set of reference genomes and computes a score for each reference genome. We used MetaNISH to scan ~300,000 SRA metagenomic runs from 2010 onwards and identified five datasets containing C. auris reads. Finally, GridRepublic has implemented a prospective C. auris molecular monitoring system using MetaNISH and volunteer computing.

Understanding the exposure risk of aerosolized Coccidioides in a Valley fever endemic metropolis.

Coccidioides is the fungal causative agent of Valley fever, a primarily pulmonary disease caused by inhalation of fungal arthroconidia, or spores. Although Coccidioides has been an established pathogen for 120 years and is responsible for hundreds of thousands of infections per year, little is known about when and where infectious Coccidioides arthroconidia are present within the ambient air in endemic regions. Long-term air sampling programs provide a means to investigate these characteristics across space and time. Here we present data from > 18 months of collections from 11 air sampling sites across the Phoenix, Arizona, metropolitan area. Overall, prevalence was highly variable across space and time with no obvious spatial or temporal correlations. Several high prevalence periods were identified at select sites, with no obvious spatial or temporal associations. Comparing these data with weather and environmental factor data, wind gusts and temperature were positively associated with Coccidioides detection, while soil moisture was negatively associated with Coccidioides detection. These results provide critical insights into the frequency and distribution of airborne arthroconidia and the associated risk of inhalation and potential disease that is present across space and time in a highly endemic locale.

A Phylogeographic Description of Histoplasma capsulatum in the United States.

Histoplasmosis is one of the most under-diagnosed and under-reported endemic mycoses in the United States. Histoplasma capsulatum is the causative agent of this disease. To date, molecular epidemiologic studies detailing the phylogeographic structure of H. capsulatum in the United States have been limited. We conducted genomic sequencing using isolates from histoplasmosis cases reported in the United States. We identified North American Clade 2 (NAm2) as the most prevalent clade in the country. Despite high intra-clade diversity, isolates from Minnesota and Michigan cases were predominately clustered by state. Future work incorporating environmental sampling and veterinary surveillance may further elucidate the molecular epidemiology of H. capsulatum in the United States and how genomic sequencing can be applied to the surveillance and outbreak investigation of histoplasmosis.

Public Health Research Priorities for Fungal Diseases: A Multidisciplinary Approach to Save Lives.

Fungal infections can cause severe disease and death and impose a substantial economic burden on healthcare systems. Public health research requires a multidisciplinary approach and is essential to help save lives and prevent disability from fungal diseases. In this manuscript, we outline the main public health research priorities for fungal diseases, including the measurement of the fungal disease burden and distribution and the need for improved diagnostics, therapeutics, and vaccines. Characterizing the public health, economic, health system, and individual burden caused by fungal diseases can provide critical insights to promote better prevention and treatment. The development and validation of fungal diagnostic tests that are rapid, accurate, and cost-effective can improve testing practices. Understanding best practices for antifungal prophylaxis can optimize prevention in at-risk populations, while research on antifungal resistance can improve patient outcomes. Investment in vaccines may eliminate certain fungal diseases or lower incidence and mortality. Public health research priorities and approaches may vary by fungal pathogen.

Application of Real-Time PCR Assays for the Diagnosis of Histoplasmosis in Human FFPE Tissues Using Three Molecular Targets.

Histoplasmosis is a fungal infection caused by the thermally dimorphic fungus Histoplasma capsulatum. This infection causes significant morbidity and mortality in people living with HIV/AIDS, especially in countries with limited resources. Currently used diagnostic tests rely on culture and serology but with some limitations. No molecular assays are commercially available and the results from different reports have been variable. We aimed to evaluate quantitative real-time PCR (qPCR) targeting three protein-coding genes of Histoplasma capsulatum (100-kDa, H and M antigens) for detection of this fungus in formalin-fixed paraffin-embedded (FFPE) samples from patients with proven histoplasmosis. The sensitivity of 100-kDa, H and M qPCR assays were 93.9%, 91% and 57%, respectively. The specificity of 100-kDa qPCR was 93% when compared against samples from patients with other mycoses and other infections, and 100% when samples from patients with non-infectious diseases were used as controls. Our findings demonstrate that real-time PCR assays targeting 100-kDa and H antigen showed the most reliable results and can be successfully used for diagnosing this mycosis when testing FFPE samples.

Comparing genomic variant identification protocols for Candida auris.

Genomic analyses are widely applied to epidemiological, population genetic and experimental studies of pathogenic fungi. A wide range of methods are employed to carry out these analyses, typically without including controls that gauge the accuracy of variant prediction. The importance of tracking outbreaks at a global scale has raised the urgency of establishing high-accuracy pipelines that generate consistent results between research groups. To evaluate currently employed methods for whole-genome variant detection and elaborate best practices for fungal pathogens, we compared how 14 independent variant calling pipelines performed across 35 Candida auris isolates from 4 distinct clades and evaluated the performance of variant calling, single-nucleotide polymorphism (SNP) counts and phylogenetic inference results. Although these pipelines used different variant callers and filtering criteria, we found high overall agreement of SNPs from each pipeline. This concordance correlated with site quality, as SNPs discovered by a few pipelines tended to show lower mapping quality scores and depth of coverage than those recovered by all pipelines. We observed that the major differences between pipelines were due to variation in read trimming strategies, SNP calling methods and parameters, and downstream filtration criteria. We calculated specificity and sensitivity for each pipeline by aligning three isolates with chromosomal level assemblies and found that the GATK-based pipelines were well balanced between these metrics. Selection of trimming methods had a greater impact on SAMtools-based pipelines than those using GATK. Phylogenetic trees inferred by each pipeline showed high consistency at the clade level, but there was more variability between isolates from a single outbreak, with pipelines that used more stringent cutoffs having lower resolution. This project generated two truth datasets useful for routine benchmarking of C. auris variant calling, a consensus VCF of genotypes discovered by 10 or more pipelines across these 35 diverse isolates and variants for 2 samples identified from whole-genome alignments. This study provides a foundation for evaluating SNP calling pipelines and developing best practices for future fungal genomic studies.

Role of Microbiota in the Skin Colonization of Candida auris.

Candida auris is an emerging multidrug-resistant fungal pathogen that can cause life-threatening infections in humans. Unlike other Candida species that colonize the gut, C. auris efficiently colonizes the skin and contaminates the patient's environment, resulting in rapid nosocomial transmission and outbreaks of systemic infections. As the largest organ of the body, the skin harbors beneficial microbiota that play a critical role to protect from invading pathogens. However, the role of skin microbiota in the colonization and pathogenesis of C. auris remains to be explored. With this perspective, we review and discuss recent insights into skin microbiota and their potential interactions with the immune system in the context of C. auris skin colonization. Understanding microbiota, C. auris, and host interactions in the skin is important to develop microbiome-based therapeutic approaches to prevent and treat this emerging fungal pathogen in humans.

Genomic description of human clinical Aspergillus fumigatus isolates, California, 2020.

Aspergillus fumigatus, an environmental mold, causes life-threatening infections. Studies on the phylogenetic structure of human clinical A. fumigatus isolates are limited. Here, we performed whole genome sequencing of 24 A. fumigatus isolates collected from 18 patients in U.S. healthcare facilities in California. Single-nucleotide polymorphism (SNP) differences between patient isolates ranged from 187 to 70 829 SNPs. For five patients with multiple isolates, we calculated the within-host diversities. Three patients had a within-host diversity that ranged from 4 to 10 SNPs and two patients ranged from 2 to 16 977 SNPs. Findings revealed highly diverse A. fumigatus strains among patients and two patterns of diversity for isolates that come from the same patient, low and extremely high diversity.

Aspergillus fumigatus is an environmental mold. It can cause a severe infection called aspergillosis in patients with weakened immune systems. We analyzed A. fumigatus DNA from patients at California hospitals. We described genetic diversity between samples from the same patients and among different patients. Our findings provide insight on using genomic sequencing to investigate aspergillosis in hospitals.

Investigation of a Candida auris outbreak in a skilled nursing facility - Virginia, United States, October 2020-June 2021.

Candida auris, an emerging multi-drug resistant organism, is an urgent public health threat. We report on a C. auris outbreak investigation at a Virginia ventilator skilled nursing facility. During October 2020-June 2021, we identified 28 cases among residents in the ventilator unit. Genomic evidence suggested ≥2 distinct C. auris introductions to the facility. We identified multiple infection and prevention control challenges, highlighting the importance of strengthening multi-drug resistant organism prevention efforts at ventilator skilled nursing facilities.

Genomic Epidemiology Linking Nonendemic Coccidioidomycosis to Travel.

Coccidioidomycosis is a fungal infection endemic to hot, arid regions of the western United States, northern Mexico, and parts of Central and South America. Sporadic cases outside these regions are likely travel-associated; alternatively, an infection could be acquired in as-yet unidentified newly endemic locales. A previous study of cases in nonendemic regions with patient self-reported travel history suggested that infections were acquired during travel to endemic regions. We sequenced 19 Coccidioides isolates from patients with known travel histories from that earlier investigation and performed phylogenetic analysis to identify the locations of potential source populations. Our results show that those isolates were phylogenetically linked to Coccidioides subpopulations naturally occurring in 1 of the reported travel locales, confirming that these cases were likely acquired during travel to endemic regions. Our findings demonstrate that genomic analysis is a useful tool for investigating travel-related coccidioidomycosis.

Genomics and metagenomics of Madurella mycetomatis, a causative agent of black grain mycetoma in Sudan.

Madurella mycetomatis is one of the main causative agents of mycetoma, a debilitating neglected tropical disease. Improved understanding of the genomic diversity of the fungal and bacterial causes of mycetoma is essential to advances in diagnosis and treatment. Here, we describe a high-quality genome assembly of M. mycetomatis and results of the whole genome sequence analysis of 26 isolates from Sudan. We demonstrate evidence of at least seven genetically diverse lineages and extreme clonality among isolates within these lineages. We also performed shotgun metagenomic analysis of DNA extracted from mycetoma grains and showed that M. mycetomatis reads were detected in all sequenced samples with the average of 11,317 reads (s.d. +/- 21,269) per sample. In addition, 10 (12%) of the 81 tested grain samples contained bacterial reads including Streptococcus sp., Staphylococcus sp. and others.

Molecular Epidemiology of Blastomyces gilchristii Clusters, Minnesota, USA.

We characterized 2 clusters of blastomycosis cases in Minnesota, USA, using whole-genome sequencing and single-nucleotide polymorphism analyses. Blastomyces gilchristii was confirmed as the cause of infection. Genomic analyses corresponded with epidemiologic findings for cases of B. gilchristii infections, demonstrating the utility of genomic methods for future blastomycosis outbreak investigations.

MycoSNP: A Portable Workflow for Performing Whole-Genome Sequencing Analysis of Candida auris.

Candida auris is an urgent public health threat characterized by high drug-resistant rates and rapid spread in healthcare settings worldwide. As part of the C. auris response, molecular surveillance has helped public health officials track the global spread and investigate local outbreaks. Here, we describe whole-genome sequencing analysis methods used for routine C. auris molecular surveillance in the United States; methods include reference selection, reference preparation, quality assessment and control of sequencing reads, read alignment, and single-nucleotide polymorphism calling and filtration. We also describe the newly developed pipeline MycoSNP, a portable workflow for performing whole-genome sequencing analysis of fungal organisms including C. auris.

Comparison between Two Molecular Techniques: Nested and Real-Time Polymerase Chain Reaction Targeting 100-kDa Hc Protein for Detection of Histoplasma capsulatum in Environmental Samples.

Histoplasmosis, one of the most frequent endemic mycoses in the Americas, is caused by the inhalation of airborne conidia of Histoplasma capsulatum. Better understanding of the distribution of this fungus in the environment is important for the development of appropriate public health measures to prevent human infections. Previously, we used Hc100 nested polymerase chain reaction (PCR) to identify H. capsulatum DNA in 10% of environmental samples in Colombia. Here, we validate a 100-kDa real-time PCR assay for the detection of this fungus in the environment. Using this method, we identified H. capsulatum DNA in 80% of samples of raw organic materials, such as chicken manure, soil from caves, and bird and bat guano, as well as in 62% of samples of organic fertilizer that underwent the composting process. We demonstrated that 100-KDa real-time PCR is a useful tool for environmental surveillance that can be used to identify the potential reservoirs of H. capsulatum and to prevent outbreaks, especially in people with the higher risk of exposure, such as spelunkers, farmers, poultry manure collectors, and anyone who handle organic fertilizers or bat and bird excreta.

Factors Influencing Distribution of Coccidioides immitis in Soil, Washington State, 2016.

Coccidioides immitis and Coccidioides posadasii are causative agents of Valley fever, a serious fungal disease endemic to regions with hot, arid climate in the United States, Mexico, and Central and South America. The environmental niche of Coccidioides spp. is not well defined, and it remains unknown whether these fungi are primarily associated with rodents or grow as saprotrophs in soil. To better understand the environmental reservoir of these pathogens, we used a systematic soil sampling approach, quantitative PCR (qPCR), culture, whole-genome sequencing, and soil chemical analysis to identify factors associated with the presence of C. immitis at a known colonization site in Washington State linked to a human case in 2010. We found that the same strain colonized an area of over 46,000 m2 and persisted in soil for over 6 years. No association with rodent burrows was observed, as C. immitis DNA was as likely to be detected inside rodent holes as it was in the surrounding soil. In addition, the presence of C. immitis DNA in soil was correlated with elevated levels of boron, calcium, magnesium, sodium, and silicon in soil leachates. We also observed differences in the microbial communities between C. immitis-positive and -negative soils. Our artificial soil inoculation experiments demonstrated that C. immitis can use soil as a sole source of nutrients. Taken together, these results suggest that soil parameters need to be considered when modeling the distribution of this fungus in the environment. IMPORTANCE Coccidioidomycosis is considered a highly endemic disease for which geographic range is likely to expand from climate change. A better understanding of the ecological niche of Coccidioides spp. is essential for generating accurate distribution maps and predicting future changes in response to the changing environment. Our study used a systematic sampling strategy, advanced molecular detection methods, and soil chemical analysis to identify environmental factors associated with the presence of C. immitis in soil. Our results demonstrate the fungus can colonize the same areas for years and is associated with chemical and microbiological soil characteristics. Our results suggest that in addition to climate parameters, soil characteristics need to be considered when building habitat distribution models for this pathogen.

Genomic Diversity of Azole-Resistant Aspergillus fumigatus in the United States.

Azole resistance in pathogenic Aspergillus fumigatus has become a global public health issue threatening the use of medical azoles. The environmentally occurring resistance mutations, TR34/L98H (TR34) and TR46/Y121F/T289A (TR46), are widespread across multiple continents and emerging in the United States. We used whole-genome single nucleotide polymorphism (SNP) analysis on 179 nationally represented clinical and environmental A. fumigatus genomes from the United States along with 18 non-U.S. genomes to evaluate the genetic diversity and foundation of the emergence of azole resistance in the United States. We demonstrated the presence of clades of A. fumigatus isolates: clade A (17%) comprised a global collection of clinical and environmental azole-resistant strains, including all strains with the TR34/L98H allele from India, The Netherlands, the United Kingdom, and the United States, and clade B (83%) consisted of isolates without this marker mainly from the United States. The TR34/L98H polymorphism was shared among azole-resistant A. fumigatus strains from India, The Netherlands, the United Kingdom, and the United States, suggesting the common origin of this resistance mechanism. Six percent of azole-resistant A. fumigatus isolates from the United States with the TR34 resistance marker had a mixture of clade A and clade B alleles, suggestive of recombination. Additionally, the presence of equal proportions of both mating types further suggests the ongoing presence of recombination. This study demonstrates the genetic background for the emergence of azole resistance in the United States, supporting a single introduction and subsequent propagation, possibly through recombination of environmentally driven resistance mutations. IMPORTANCE Aspergillus fumigatus is one of the most common causes of invasive mold infections in patients with immune deficiencies and has also been reported in patients with severe influenza and severe acute respiratory syndrome coronavirus 2 (SARs-CoV-2). Triazole drugs are the first line of therapy for this infection; however, their efficacy has been compromised by the emergence of azole resistance in A. fumigatus, which was proposed to be selected for by exposure to azole fungicides in the environment [P. E. Verweij, E. Snelders, G. H. J. Kema, E. Mellado, et al., Lancet Infect Dis 9:789-795, 2009, https://doi.org/10.1016/S1473-3099(09)70265-8]. Isolates with environmentally driven resistance mutations, TR34/L98H (TR34) and TR46/Y121F/T289A (TR46), have been reported worldwide. Here, we used genomic analysis of a large sample of resistant and susceptible A. fumigatus isolates to demonstrate a single introduction of TR34 in the United States and suggest its ability to spread into the susceptible population is through recombination between resistant and susceptible isolates.